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Abstract 

It is shown that the recently established strong mass-dependence 
of the radii of the hadron sources, as observed in HBT analyses of the 
e + e~ annihilation, can be explained by assuming a generalized inside- 
outside cascade, i.e. that (i) the /crar-momenta and the space-time 
position four-vectors of the produced particles are approximately pro- 
portional to each other and (ii) the "freeze-out" times are distributed 
along the hyperbola t 2 — z 2 = Tq . 

It has been found recently that the parameters describing the B-E inter- 
ference in e + e~ annihilation depend strongly on the masses of the particles 
used in the analysis ^ |2|. One finds r n between 0.7 and 1 fm; between 
0.5 and 0.7 fm; between 0.1 and 0.2 fm. 

In the present note we suggest that this dependence can be understood 
if the produced particles satisfy approximately the (generalized) Bjorken- 
Gottfried conditions fBL ffl : 
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(i) The 4-momentum and the 4-vector x M describing the space-time 
position of the production ("freeze-out") point of a particle are proportional 



(1) 



The proportionality factor A is a scalar with respect to boost in the longitu- 
dinal direction; 

(ii) Particles are produced at a fixed proper time r after the collision, 

t 2 ~z 2 = rl (2) 

where t, z are time and longitudinal position of the production point. 
From (|l|) and (0) we derive 



A 



M ± 

To 



(3) 



where M\ = E 2 — gj?. Thus finally we have 



To 



-x. 



(4) 



This picture is, of course, purely classical and can only be treated as a 
heuristic guide-line when applied to actual production processes. A more 
adequate formulation of these conditions can be achieved using the Wigner 
representation W(P, x) of the (single-particle) density matrix which, as is well 
known (see e.g. 0), corresponds - as close as possible without contradicting 
quantum mechanics - to the space-time and momentum distribution of the 
produced particles. To implement the conditions (i), (ii) above, we postulate 
W(P,x) in the form 



W(P, x) ~ 5(t 2 -z 2 - t 2 ) exp 



.1 I 



2R 2 



PL 
2A 2 ± 



exp 



(P + - ^* + ) 2 + (P. - ^x_) 2 (Px-^ 



where 
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x± — t ± z\ 
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so that 

M\ = P + P_; r 2 = (7) 

The first exponential represents a standard cylindrically symmetric "longitu- 
dinal" distribution in momentum and in configuration spaced. The new point 
is the second exponential which introduces correlation between the momen- 
tum and the point of emission of the particle, as required by the generalized 
Bj or ken-Gottfried condition (Q). Such correlations are known to influence 
strongly the HBT effect on particle spectra ||. It is thus this factor! which, 
we think, is responsible for the mass dependence of the observed HBT radii. 

To derive HBT correlations we need to calculate from (|5]) the density 
matrix in momentum space (see e.g. |8|, [| [KJ). This can be done using the 
relation between W(P, x) and p(q, q') which reads 



P- 



Q 

2 ' 



J d^xe iQx W{P, x) 



(8) 



From (m we see that now we have to take 



P 



q + q' 



Ml 



Q 



(9) 



with the 4-momenta q and q' on the mass-shell. 

To continue, it is convenient, as usual, to introduce the rapidities 



n 



1 X A 

2 x. 



(10) 



The longitudinal integral 



dijexp 



M 



252 



exp (itq [m± cosh(y — rj) — m' ± cosh(y' — r/)}) 



(11) 



lr To simplify the argument, we ignore the longitudinal momentum and z dependence 
of the single particle spectrum. This seems a reasonable approximation at high energy. 

2 Admittedly, the form (^) is rather schematic. In particular, gausssians are taken for 
simplicity and can be replaced if necessary. We also did not include fluctuations of To. 
These simplifications are not essential for our argument, however. 
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where (m±,y) and (m' ± ,y') are transverse masses and rapidities correspond- 
ing to momenta q and q', can be approximated by 



dr/exp 



Ml 



exp itq 



m± 



\ (v-yf 



m 



1 (v - y')' 



(12) 



Ignoring normalization and phase factors, inessential for our argument, we 
thus obtain 



/ M 2 

exp ( -ir ^(m ± (Y - yf - m' ± (Y - y'f) 



exp (^-^m±m ± (y - y'f^j (13) 



where 
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Ml T 

— ^ i — (m± — m I 



(14) 



The transverse integral can be evaluated exactly. Ignoring again the normal- 
ization and phase factors we have 



I± = exp 



P 2 



/ 



2A : 



d 2 x exp 
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(15) 



where all vectors are two-dimensional (transverse) and 



A 2 



+ 



R 2 ± r 2 5 2 ± 



MlR\ + r 2 5 2 ' A 2 ' 



p2 _ 



(16) 



From (JT5I) we find the single particle transverse momentum distribution: 



da 
d 2 q± 



I.l(q= $ = 9l) = exp 



o 2 r 2 



'2(m 2 + g 2 ) J R 2 L + 2r 2 (5 2 2A^ 



(17) 
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One sees that the average transverse momentum is largely determined by the 
value of A± which thus cannot be too large if one wants to insure average 
transverse momentum smaller than, say, 500 MeV. 

Let us also note at this point that consistency with uncertainty principle 
implies the inequality JTOf 

RefAeff > I (18) 



As seen from (jig) , at large transverse mass M±, this inequality can only 
be satisfied if 8± is significantly larger than Aj_ (and thus than the average 
transverse momentum). 

To proceed, we shall assume that all correlations between particles which 
are not caused by Bose-Einstein interference can be neglected. Using the 
formulation of (IIJ we thus write the two-particle density matrix as a product 

p(qi, Q2\ q[, q' 2 ) = p(?i> i[)p(q2i q' 2 ) (!9) 

It then follows from the general theory of HBT effect (see, e.g. |§) that the 
observed two-particle distribution is given by 

Q2) = p(qi,ii)p(<i2, g 2 ) + p(qi, q 2 )p{q 2 , gi) = ^(gi)^(g 2 )(l ± C(q u g 2 )) 

(20) 

where 



C(gi,g2) = C|i(gi,g2)c±(gi,g 2 ) 
\h(qi,q2)\ 2 \i±(qi,q2)\ 2 
h gi)^ii(g2, g 2 ) gi)i±(g 2 , g 2 ) 



(21) 



describes the HBT correlations. 
Using (p~5|) we find 

C ± = e -(^-fc) 2R HBT = e -Q\ R \ H BT (22) 

where 

It 2 11 

RIhbt = Riff ~ = ^x^^/i - 4) - TaT ( 23 ) 



Since 



Ml = ( m ^ + m ^ ) 2 + mi±m2± sinh2 (yi^Vl) (24 ) 
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we conclude that indeed R\hbt ^ a ^ s with increasing (transverse) mass of the 
particle. 

For Cm we have 



C||(<?i)92) = exp [-fl|| MT (mu|/i - m 2 ±y 2 + (m lx - m 2 ±)Y) j (25) 
where 

r 2 M 2 

R " HBT = 2\B^I (26) 

From (0), one sees that also R^ HBT falls with increasing M\. 

If, as is customary (see e.g. JT|), one works in the frame where Y — 0, 



5f) can be written as 

C\\(qi, q 2 ) « e " fl ii« OT(<?lz " 92z)2 = e^** 3,0 " (27) 

This completes the qualitative discussion of the mass effect in our ap- 
proach. It remains to be seen if the values of the HBT radii given by (^) 
and fl26|) can be adjusted to be close to the ones obtained from the LEP data 



In Fig. 1 R\\hbt an d R±hbt are plotted versus M±, the transverse mass 
of the two-particle system. The values of other parameters were taken as 
follows : A ± = 360 MeV, r = R ± = 1.2 fin, 5 ± = 700 MeV, S\\ = 350 MeV, 
[mix — Tn,2±\ = 150 MeV. One sees a rather strong mass dependence of both 
longitudinal and transverse radii. We did not try to fit the obtained values to 
the data as this would require working directly with data themselves and thus 
goes beyond the scope of the present investigation. It is nevertheless recom- 
forting to observe that the HBT radii, obtained with "reasonable" values of 
the model parameters, are not far from the ones found in LEP experiments. 

We thus conclude that the existing data on HBT radii are consistent with 
the hypothesis that -in e + e~ annihilation at high energy- 4-momentum of a 
produced particle is approximately proportional to its space time position 
4-vector at the freeze-out time!. 

This proportionality is of course well-known for the longitudinal compo- 
nents [0, |J, and is exhibited explicitly in numerous models ||TT|. At this point 



our approach is similar to the one proposed for a longitudinally expanding 



^Recently an alternative interpretation has been proposed in Ji2[ 
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T„ = R A - 1.2 fm 
R HB T tfml a 1= 360 MeV 
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Figure 1. R±hbt and R\\hbt plottet versus M±. The parameters are shown 
in the figure. The data from nir, KK and AA correlations are also indicated. 



fireball [[L3| |9|, although the mass dependence following from our Eq. (jg) 
seems somewhat stronger. On the other hand, a rather novel feature follow- 
ing from our analysis is that the original Gottfried- Bj or ken proportionality 
relation should be extended to include also the transverse components of the 
4- vectors, as explicitely expressed in (H). 
Several comments are in order. 

(i) It should be emphasized that our argument is only semi-quantitative 
and can be improved in many details when applied to real data. In par- 
ticular, the gaussians in the Wigner function (gj) can be replaced by more 
realistic functions for numerical analysis. Also, the Fourier transform (O) 



can be calculated numerically without approximations shown in (|i2|) , which 
were introduced simply to obtain an analytic result. Finally, including a dis- 
tribution of t is probably needed to obtain a good description of data. We 
feel, however, that all this necessary fine tuning does not invalidate our main 
conclusion, summarized in Eq. (|4]). 

(ii) As we already mentioned, the results shown in Fig. 1 do not represent 
a fit to experimental data which we think would be premature at the present 
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stage. Therefore, the values of the parameters used to produce this figure are 
by no means final. Some of them seem rather stable, however. In particlular, 
Aj_ is closely related to the average transverse momentum and thus cannot be 
arbitrarily changed. Also a rather large value of 5± seems necessary to satisfy 
the consistency condition fll8|). This means that the correlation between the 
tranverse momentum and transverse position of a particle at freeze-out is 
fairly weak. It is remarkable that such a weak correlation is sufficient to create 
a strong variation of R±hbt with the transverse mass of the investigated two- 
particle system. 

(iii) From the point of view of data analysis, our argument emphasises 
the importance of the investigation of the HBT correlations as function of 
the transverse mass of the pion pair. 

(iv) Relation (|j), when applied to transverse directions, implies the exis- 
tence of an important "collective transverse flow" in the system of particles 
produced in e + e~ annihilation^!. It would be interesting to search for other 
evidence of such a "flow" in the data. 

(v) A natural modification of the relation (0) is to consider freeze-out 
times given by the fully Lorentz-invariant formula 

e-z*-y*-x* = T* (28) 

which leads to qualitatively similar results as those discussed in the present 
paper. It is not clear if the present data can distinguish between (|2]) and ( p8|) 
but investigation of this question is certainly a challenging issue for future 
work. 



vi) The recent data of L3 coll. [15] show a strong dependence of the trans- 



verse 7T7T HBT radius (and a somewhat weaker dependence of the longitudinal 
radius) on the average transverse mass of the two pions m± = |(mi_|_ + m 2 A_). 
This seems not inconsistent with our results, although more work is needed 
to establish a closer connection between M± and the average trasverse mass 
m_i_ which is used to parametrize the data. Thus before more detailed cal- 
culations (including a realistic single particle distribution) are performed, 
it is not clear to what extent the results shown in Fig.l are related to the 



observations of |15 
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